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Sonic booms and diffusion wakes generated by a heavy quark in thermal AdS/CFT 
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We evaluate the Poynting vector generated by a heavy quark moving through a thermal state of 
M = 4 gauge theory using AdS/CFT. A significant diffusion wake is observed as well as a Mach 
cone. We discuss the ratio of the energy going into sound modes to the energy coming in from the 
wake. 
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I. INTRODUCTION 

The phenomenon of jet-splitting observed at the Rel- 
ativistic Heavy Ion Collider (RHIC) 0, 0] motivates a 
study of the stress tensor generated by a quark mov- 
ing through a thermal state of Af = 4 super- Yang-Mills 
theory (SYM). Although SYM is significantly differ- 
ent from quantum chromodynamics (QCD), it is widely 
hoped that comparisons between the two will lead to in- 
sights not readily extracted from conventional techniques 
such as perturbation theory and lattice simulations. The 
advantage of studying SYM is that powerful string the- 
ory methods exist to calculate gauge-singlet observables, 
such as the stress tensor, at strong coupling, starting 
from gravitational calculations in five-dimensional anti- 
de Sitter space (AdSs). These methods rely on the fact 
that SYM is a conformal field theory (CFT). 

The basic stringtheory setup is well reviewed in earlier 
works B I, i, UU 1 3- Briefly, an infinitely massive, 
fundamentally charged quark that is dragged at constant 
velocity through an infinite, static, thermal plasma of 
SYM can be thought of as having a string trailing be- 
hind it down into AdSs-Schwarzschild, generating a drag 
force as computed in [l(| UH (see also [H, HH). This 
string perturbs the metric g^ v of AdSs-Schwarzschild, 
and those perturbations can be translated into the ex- 
pectation value (T mn ) of the stress tensor in the bound- 
ary gauge theory. Although it would perhaps be prefer- 
able to study an energetic gluon propagating through 
the thermal plasma, the infinitely heavy quark case is a 
textbook exercise in AdS/CFT because the quark has no 
dynamics apart from its constant-velocity motion: thus 
its trajectory can be regarded as part of the boundary 
conditions specified in the gravitational calculation. In- 
falling boundary conditions are specified at the horizon 
of AdSs-Schwarzschild, corresponding to causal physics 
in the boundary gauge theory. 



In the first work on jet-splitting in AdS/CFT 0], 
Fourier coefficients (T mn (k)) were computed numerically, 
and high-fc and low-fc asymptotics were given. Subse- 
quent work [EH, 0,9 has considerably refined the high-fc 
estimates. The position-space stress tensor can in prin- 
ciple be recovered from its Fourier coefficients using 



(T mn (t,x)) 



d 3 k 
(2^)3 



iki (x 1 — vt)-\-ik2X 2 -\-ik3x 3 frpmn 



Here and below, 
the positive x 1 
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(T mn (k)} . 

v is the quark's velocity directed along 
direction. Evidence of a sonic boom 
was already exhibited in [4], and in [8|, [g[ it was made 
more concrete by evaluating the full position-space en- 
ergy density (T 00 (t, x)). A match to linearized hydrody- 
namics with 77/s = 1/47T was remarked upon in [S|, |9f| - 
What AdS/CFT provides beyond hydrodynamics is an 
all-scales description of dissipation, from the near-field of 
the quark all the way out to the linearized hydro regime. 

A crucial question for application to RHIC phe- 
nomenology is the relative strength of the sonic boom 
and the diffusion wake. In the lab frame, the diffusion 
wake is a structure primarily visible in the Poynting vec- 
tor, and it is concentrated within a parabolic region be- 
hind the quark. It is closely related to the well-known 
laminar wake of non-relativistic hydrodynamics. The 
forward flow in the diffusion wake should lead to en- 
hanced particle production approximately collinear with 
the hard parton, as compared to vacuum fragmentation. 
The sonic boom, in contrast, describes high-angle emis- 
sion. Intriguingly, data from the PHENIX experiment 
PI favors suppression of the diffusion wake to such an 
extent that earlier authors [l4| have suggested turning 
it off altogether in a linearized hydro analysis. Other 
authors have proposed a model 15] in which it is as- 
sumed that 75% percent of energy lost goes into a sonic 
boom. Data from STAR does not show as clear a min- 
imum in associated particle production in the direction 
of the away-side parton [2|. In the string theory calcu- 
lation, the relative strength of the sonic boom and the 
diffusion wake is fixed. Because other first-principles es- 
timates of the relative strength of these two effects are 
lacking, it is clearly of interest to work out the string the- 



2 



ory predictions in some detail. All the pieces are already 
in place except for sufficiently precise high-fc estimates 
of (T° 4 (fc)) and a robust numerical scheme for perform- 
ing the Fourier transforms. In this letter, we fill in these 
last gaps and extract explicit results for the Poynting 
vector. We also discuss one way of quantifying the rela- 
tive strength of sonic and diffusive modes which suggests 
that, in the rest frame of the plasma, the sonic boom 
dissipates I + v 2 times as much energy as the diffusion 
wake feeds in. 



II. DEFINITIONS 

The only scale in the problem is the temperature, so 
we will consistently work with dimcnsionlcss quantities 



X 
K 



k/irT . 



(2) 
(3) 



Instead of computing the expectation value of the com- 
plete stress tensor, we focus on a rescaled, subtracted 
quantity 



r mn {x) 



{■*T)^g 2 YM N 

/ r pmn\ — 7T 2 ( at 2 
[-L /bath = -5-(iV 



(V»"(0,z)) - (T ro ") bat h)(4) 
-l)T 4 diag{3, 1,1,1}. (5) 



On the right hand side of ^ we have set t = 0, which 
leads to no loss of information because the computa- 
tions are all done in a steady state approximation. The 
powers of irT render T mn dimensionless. The factor of 
1 / yJg'Y M N cancels out an overall scaling with \Jgy M N 
of all disturbances due to the trailing string. This scal- 
ing arises because it's how the string tension depends on 
gyu and N. The quantity (T m ")b a th is just the stress 
tensor of the SYM plasma in the absence of the quark. 
Our focus is on the first row of T" m , namely 



T 



Om 



(£ Si 6>2 S3 



(6) 



where £ and S are the rescaled, bath-subtracted energy 
density and Poynting vector. (We work in signature 
— I — I — (-, so there is no distinction between S l and Si.) 
We further define 



E 
S 



£ ^Coulomb 
S — tScoulomb ■ 



(7) 
(8) 



where fcouiomb and <Scouiomb characterize the Coulombic 
near-field of the quark. Thus E and S exclude not only 
the contribution of the bath in the absence of the quark, 
but also the contribution of the quark in the absence of 
the bath. The quantities which we eventually plot are E, 
Si , and S± as functions of Xi and X± . Here 



X± = 



(9) 



III. ASYMPTOTICS AND SUBTRACTIONS 

The key to a robust numerical evaluation of stress- 
tensor components in position space is to make a split 



£ = £ vw 
S = &™ 

cUV 



£ m - 



£rcs 

- s rcs 



(10) 

(11) 



where £ uv and S vy are analytical approximations to 
£ and S for large K which are valid to order K~ 3 , 
£ IR and S 1 ^ are analytical approximations for small K 
valid to order K°, and the remainders £ ICS and S rcs are 
uniformly bounded and integrable, so that they can be 
passed through a fast Fourier transform (FFT) with con- 
trollable errors. The full expressions for the analytical 
approximations are too long to be presented in full, but 
the leading terms are 
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48^K'f{l-v 2 )- 




Muv 


5 1 uv = -| v / ^ 2 (l-« 2 )- 
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2vKf(l-v 2 ) 


+ 3i>Muv 




-Kl + 




5 uv _ vKiK 2 (l - 
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OiR- 1 ) (12) 



£ m = - 



24 v /i?F(l^ 
1 ZivK x (l 



V 



K 2 ± 
^ - 3v 2 K 2 



+ 0{K- 1 ) 
^OiK- 1 ) 



Mir 



IR 



Si 



Sl R 



2tt K 2 - 3v 2 K 2 - ivK 2 Ki 
1 ZivKijl + v 2 ) - 3v 2 K 2 
2^K 2 - 3v 2 K 2 - ivK 2 Ki 

1 i(l + v 2 )Ki + vK 2 - 2v 3 K 2 
~2^ K 2 - 3v 2 K 2 - ivK 2 Ki 
2v 1 + iKx/Av 



7T K 2 - 

1 i(l 



AivKi 
^v 2 )K 



+ (regulators) 
2 + v 3 KiK 2 



2tt K 2 - 3v 2 K 2 - ivK 2 Ki 
1 iK 2 



2ir K 2 - UvKi 



+ (regulators) . 



(13) 
(14) 

(15) 
(16) 
(17) 



and similarly for K± and Sj 



Expressions for S^ v and iS| R may be deduced by rota- 
tional invariance around the K\ axis. The omitted terms 
in (fll)|) and (JT7J) denoted "(regulators)" are analogous to 
the second term in (fl"5)) . 

The strategy for obtaining UV subtractions was ex- 
plained in [8(: we solve for the metric perturbations at 
large K using an iterative Green's function approxima- 
tion, as first introduced in 0, then "soften" the result- 
ing power series through the introduction of a dimen- 
sionless parameter /xyv as in (|12p - |fL4"|) . To obtain the 
IR subtractions, we draw upon results of [1, [|| to find 
series expansions at small K for components of 7" m ™, 
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then partially resum the series, as in (|15 p - (fT7|) , to elim- 
inate divergences on the Mach cone and the Xi axis, 
then add Pauli-Villars-style regulators involving another 
dimensionless parameter ^ir. 

The terms in (II 5[) — (| 1 T[) with cubic denominators do not 
appear to admit analytical Fourier transforms; however 
they can easily be Fourier transformed in the K± di- 
rections, leaving a final one-dimensional FFT to be per- 
formed numerically. Some of the subleading terms in 
our subtractions do not by themselves satisfy appropri- 
ate conservation properties, so the conservation of the full 
stress tensor has to be restored by combining numerical 
and analytical results. 



and 1 + v 2 times the energy fed in by the diffusion wake 
[lij . These ratios should be regarded with some caution 
because the argument leading to them is essentially al- 
gebraic; however, integrating the energy flux across the 
diffusion wake confirms (|24|) . 

Instead of calculating in the rest frame of the plasma, 
as we have done, one may instead perform computations 
analogous to (f2"3"]l and in the rest frame of the quark. 
Then, using also the IR asymptotics of T%% obtained in 
[H, one can show that the energy lost by the quark into 
sound modes precisely cancels the energy fed in through 
diffusion modes. Some intuition regarding this result may 
be gained by noting that in the quark rest frame, the zero 
component of the drag force vanishes. 



IV. CONSERVATION 



The stress tensor must be conserved except at the lo- 
cation of the quark, and the failure of conservation there 
matches the drag force on the quark [i[ . In terms of our 
rescaled variables, and in the rest frame of the plasma, 







-r mn = j n = -rg &s 8{x l - vx°) 



dX" 



After using the co-moving ansatz, (fT9|) becomes 
K m = (-vKx K x K 2 K 3 ) . 



(18) 
(19) 



(20) 
(21) 



Now consider a split of the IR asymptotic values of T mn 

(22) 



/rOm r-r-Qm 

IR sound 



diffuse 



where we define T®ound as a ^ the terms in (fI5|l - (fI7|> with 
cnbic polynomials in the denominator; ^°iff usc by defini- 
tion comprises all the remaining terms. With this split, 



lim iK m T?r 



v 



'■m- 4 diffuse 



2tt 
1 

'2tt 



= -(1 + ^)^(23) 



>•> dr; 



(24) 



The m = term in the left hand side of (|23|) is, for- 
mally, the time derivative of the total energy in sound 
modes. Adding in the m = 1, 2, 3 terms accounts for en- 
ergy flow into sound modes at large X, and the sum over 
all m is the total rate of energy loss by the quark into 
sound modes. Likewise, (|24|) is the energy loss into diffu- 
sion modes. But the relative sign between (|23p and (j2~4")l 
shows that while sound modes carry energy away from 
the quark, the diffusion wake actually feeds energy in 
toward the quark. 

Energy loss to sound modes, as measured by (|2"3"|) , is 
greater by a factor of 1 + l/v 2 than the total energy loss, 



V. NUMERICAL RESULTS 

Our numerical analysis was mostly done using 128 3 
grids for the three-dimensional FFT's of £ rcs and <S ros , 
with wave-numbers K ranging from —10 to 10 in each 
component. The one-dimensional FFT's of the terms in 
£ IR and S 5 ^ with cubic denominators were performed on 
finer grids, usually with 1944 points, with K\ ranging 
from —20 to 20. We chose /^uv = Wr = 1- Conser- 
vation was checked by comparing £ to V • S in position 
space. This is a stringent check because derivatives of 
numerically known functions tend to be noisy. Conser- 
vation worked to within a few percent of £ except when 
£ became very small. A summary of our numerical re- 
sults for v = 0.75 is shown in figure [T] A more complete 
presentation of our numerical results is available [17] . 

VI. CONCLUSIONS 

The Poynting vector describing the disturbances of a 
thermal plasma of M = 4 super- Yang-Mills theory from 
a moving quark exhibits a sonic boom and a diffusion 
wake, as expected on general grounds. Using a natural 
split of the energy flow into sonic and diffusive modes, 
one finds — in the rest frame of the plasma — that energy 
is lost from the quark through sound modes and fed in 
toward the quark through the diffusion wake. These two 
effects stand in the ratio 1 + v 2 : —1. Thus the sonic 
boom is a strong effect, but the diffusion wake is also 
significant. 
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FIG. 1: Energy density and components of the Poynting vector. Sp means S± = \/ S% + S| . \X\ = 1 corresponds to 0.2 fm 
from the quark if T ~ 318 MeV. Red means positive, blue means negative, and white is close to zero. The lower right plot is 
colored according to the magnitude of S, while the arrows show the direction of S. For Xi > 3, S is so small that we probably 
cannot determine its direction reliably. The Mach cone is shown in dashed green, and the solid blue line shows where the large 
distance analytic estimate of the diffusion wake profile reaches half its maximum value. 
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